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ABSTRACT 

We report results of VRI photometric monitoring of the T Tauri star plus disk system HH30. We find 
that HH30 is highly variable over timescales of a few days with AV^ ~ 1.5 mag, A/ ~ 1.1 mag. Further- 
more we find hints of periodicity with periodograms indicating possible periods of 11.6 and 19.8 days. 
The VRI photometry is available through the anonyinous ftp service. We model the variability with 
Monte Carlo radiation transfer simulations for a spotted star plus disk system and find that the large 
variability is best reproduced with a single hot spot and circumstellar grains that are larger than typical 
interstellar grains. The apparent existence of a single hot spot and the need for larger, more forward 
throwing grains is consistent with previous modeling of HST imagery. 

Subject headings: stars: pre-main-sequence — stars: rotation — stars: spots — stars: individual 
(HH30) — accretion, accretion disks — radiative transfer 



1. INTRODUCTION 

Signatures of hot and cool starspots have been observed 
in numerous T Tauri stars, primarily in the form of pho- 
tometric variability, with periodic variability detected in 
hundreds of systems (e.g., Bouvier et al. 1993; Wichmann 
et al. 1998; Choi & Herbst 1996; Makidon et al. 1997; 
Stassun et al. 1999). While cool spots have been mapped 
through Doppler imaging of weak- lined T Tauri stars (e.g., 
Hatzes 1995), hot spots tend to be associated with vari- 
ability in Classical T Tauri stars (e.g., Kenyon et al. 1994; 
Bouvier et al. 1993; Herbst et al. 1994). The hot spots are 
thought to be associated with the accretion process, and 
the currently favored magnetic accretion models naturally 
provide for hot spots on stellar surfaces (Ghosh & Lamb 
1979a, 1979b; Konigl 1991; Shu et al. 1994; Ostriker & 
Shu 1995; Najita 1995). 

Hot spots yield specific brightness and polarization vari- 
ations as they rotate into and out of view (Wood et al. 
1996; Mahdavi & Kenyon 1998; Stassun & Wood 1999). 
Wood & Whitney (1998) investigated the effects of hot 
starspots on the morphology of scattered light disks. The 
non-axisymmetric illumination of the disk by hot spots 
leads to an asymmetric brightening of the disk. Such 
asymmetric brightening has been detected in HST obser- 
vations of the edge-on disk system HH30 (Burrows et al. 
1996; Stapelfeldt et al. 1999a), hinting at possible agree- 
ment with the magnetic accretion model. 

The gross consistency between the existing observa- 
tions of HH30 and the predictions of the magnetic accre- 
tion model is tantalizing, but ambiguous due to the large 
timescales separating the HST observations. Periodic vari- 
ability due to magnetic accretion is expected to match the 
rotational periods for Classical T Tauri stars: of the order 
less than one day to tens of days (Attridge & Herbst 1992; 



Bouvier et al. 1993; Edwards et al. 1993; Eaton et al. 
1995; Stassun et al. 1999; Herbst, et al. 2000). In some 
cases the hotspots are stable over many rotation periods. 

Motivated by the apparent success of the magnetic ac- 
cretion model at explaining the morphological variations 
observed in the limited HST imagery, we have under- 
taken a ground-based photometric monitoring campaign 
on HH30 to determine the timescale of any variability. In 
Section 2 we present the results of our observational cam- 
paign, in Section 3 we present models for the observed 
photometric variability, and we summarize our findings in 
Section 4. 

2. OBSERVATIONS 

Data for HH30 (RA = 04 31 37.5, DEC = 18 12 26.0; 
V ~ 19.5, / ~ 17.3, Mundt & Fried 1983) were ob- 
tained with the 1.2 m telescope at the F. L. Whipple 
Observatory (FLWO), using the "4Shooter" CCD mosaic 
(Szentgyorgyi et al. 2000) with four thinned, back-side 
illuminated, AR coated Loral 2048^ CCDs. The cam- 
era has a pixel scale of 0.335 arcsec/pixel and field of 
view of roughly 11.5 arcmin for each chip. The data 
were taken in the 2x2 CCD binning mode. To obtain 
multi-wavelength variability information we used VRI fil- 
ters from the FLWO "Harris Set" . We obtained data (see 
Figure 1), between September 7th 1999 and February 28th 
2000, usually consisting of up to three sets of 300 sec, 
180 sec, 180 sec VRI exposures. In total we have 24 use- 
ful images of HH30 in the F-band and 27 useful images 
in each of the i?/-bands, spanning about 174 days. The 
seeing ranged from FWHM = 1.3 arcsec to about 4.0 arc- 
se c. The VRI data can be accessed from anonymous ftp 
at [ftp : //cfa-ftp .harvard. edu/pub/kstanek/HH30 . 

Preliminary processing of the CCD frames was done 
with the standard routines in the IRAF-CCDPROC pack- 
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age.® Photometric variability was determined using two 
separate techniques: aperture photometry and point- 
spread function (PSF) fitting. PSF fitting was done 
using the photometric pipeline of the project DIRECT 
(Kaluzny et al. 1998; Stanek et al. 1998), based on 
the Daophot/ Allstar photometry package (Stetson 1987, 
1991). The differential light curves obtained by both meth- 
ods agreed well, which was reassuring considering that the 
HH30 nebulosity is somewhat resolved in our images. As 
a further check, we saw no correlation between the derived 
photometric changes and the seeing or sky brightness. 

The resulting VRI differential light curves are shown 
in Figure 1 and are adjusted so the faintest measurement 
in each band corresponds to Amag= 0. The HH30 nebu- 
losity was observed to vary in brightness over timescales 
of a few days by Ay ~ 1.4 mag, A.R ~ 1.0 mag, and 
A7 ^ 1.1 mag. The changes in brightness are very well 
correlated between the three photometric bands. There 
are hints of periodicity with the rise and fall at the start 
of our observations and again around HJD 1570 in Fig- 
ure 1. Independently for each band we searched for 
periodicity using a variant of the Lafler-Kinman (1965) 
string-length technique, proposed by Stetson (1996) and 
also the Lomb (1976) normalized periodogram method. 
The strongest signals in our periodogram techniques were 
found at 19.8 days (97% probability) and 11.6 days (95% 
probability). An 11.6 day period is more typical of T Tauri 
stars and further monitoring of HH30 will possibly yield a 
more reliable period. We emphasize that this period de- 
tection is very tentative given the present data, and we 
encourage follow-up observations to confirm or refute our 
finding. 

Photometric variability in excess of 2 mags may be ob- 
tained by hot spots on the stellar surface (e.g., Mahdavi 
& Kenyon 1998). Hot spots produce larger photometric 
variability at shorter wavelengths due to the increased 
spot/star luminosity ratio. Our observations show that 
ISV > Al, but AR < Al. The apparent inconsistency 
of the R band variability may be attributed to contami- 
nation from strong Ha and [S ii] emission associated with 
HH30's jet (Burrows et al. 1996). Bouvier et al. (1999) 
discounted a hotspot model for AA Tau because they did 
not detect any color variations (A mag^ 1.4 in all BVRI 
bands). Our HH30 data do show color variations and in 
the next section we investigate hot spot models for HH30 
and restrict our simulations to the V and / bands due to 
the likely contamination of the R band observations by jet 
emission. 

3. MODELS 

The large variability we observe may occur for a model 
with two diametrically opposed hot spots when the second 
spot is occulted (either by the star or the circumstellar 
disk) and the observed spot is seen directly at maximum 
light and is obscured by the star at minimum light. How- 
ever, for the high inclination of HH30 {i > 80°, Burrows et 
al. 1996) we only see a scattered light nebula and do not 
see the star directly. For a model with two diametrically 
opposed hot spots, we found that at high inclinations the 
photometric variability was only around 0.2 mag (Wood & 
Whitney 1998). As one spot rotates behind the star, the 
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second spot rotates to the front yielding the small ampli- 
tude variability. In addition, because we are only detecting 
scattered light we always see some reflected light from the 
hot spots yielding small variations in the total intensity. 

The fact that we are viewing scattered light makes it im- 
possible for us to apply inversion techniques to determine 
spot parameters (e.g., Bouvier et al. 1993; Vrba, Herbst, 
& Booth 1988), since the light curves are sensitive not only 
to the star and spot parameters, but also the disk struc- 
ture and circumstellar dust properties. The new WFPC2 
observations presented by Stapelfeldt et al. (1999a) and 
the large variability present in Figure 1 suggests that one 
hot spot dominates the nebular and photometric variabil- 
ity. We now construct single hot spot models for HH30. 

HH30 has been modeled by Burrows et al. (1996) and 
Wood et al. (1998) with a flared disk geometry, 

p = poexp~[z/h{zu)]'^/w" , (1) 

where w is the radial coordinate in the disk midplane and 

the scale height increases with radius, h = ho {w/R^,)^ . 
Following Burrows et al. (1996) and Wood et al. (1998) 
we adopt: /3 = 9/8, a = 15/8, ho — 0.05i?* (giving 
/i[100AU] = 15AU), and a disk mass of 2.5 x lO-^'M©. We 
construct scattered light models using our Monte Carlo ra- 
diation transfer code (Whitney & Hartmann 1992; Wood 
& Whitney 1998; Stassun & Wood 1999) and assume there 
is one circular hot spot of radius 9s at latitude (j)s on the 
stellar surface. The relative number of photons released 
from the spot and star is, 

TV, ^ 1-cosg, BxjTs) 

i + cose^Bxin) ' ^' 

where we assume that the luminosities of the star and spot 
are Planck functions at the spot (Tg) and stellar (T*) tem- 
peratures. 

Figure 2 shows the simulated V and / band photometric 
variabihty for ds = 20°, 0^ = 65°, Ts = lO'^K, T„ = 3800K 
(Wood & Whitney 1998; Stapelfeldt et al. 1999). We 
assumed that the circumstellar dust can be characterized 
by a Kim, Martin, & Hendry (1994, hereafter KMH) size 
distribution, typical of grains in the interstellar medium. 
For KMH grains the V and / band dust parameters (opac- 
ity, albedo, peak polarization, and asymmetry parameter 
in the Heyney-Greenstein scattering phase function) are 
displayed in Table 1. 

Viewing this model at low inclinations we can easily 
obtain AV ~ 1.5 mag. However at i = 82°, we cannot 
reproduce the observed large amplitude variability — we 
only sec scattered light and always detect reflected light 
from the hot spot, irrespective of rotational phase. There- 
fore the amplitude of variability is smaller (AV" ^ 0.5 mag, 
Al ~ 0.3 mag) than when we see the star directly. One 
way of increasing the variability is to change the spot size 
and temperature. Alternatively, if the forward beaming 
of the scattered light is increased (increasing the asym- 
metry parameter, </), we will see less scattered light from 
the spot at minimum light when it is occulted by the star. 
Such an increase in the asymmetry parameter, indicative 
of large grains in the disk, was found by Burrows et al. 
(1996). Other modeling efforts have also found evidence 
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for larger, more forward throwing grains in protostellar 
environments (e.g., Lucas & Roche 1998). We have there- 
fore constructed models with larger, more forward throw- 
ing grains (see parameters in Table 1). This population 
of grains was used by Cotera et al. (2000) and Whitney 
& Wolff (2000) to model multiwavclongth HST images of 
HH30. The dust grains are homogeneous spheres com- 
posed of either amorphous carbon (BEl, Roideau & Mar- 
tin 1991) or revised astronomical silicate (Wcingartner & 
Draine 2000) . The size distribution for each component is 
specified using a power-law with exponential decay (i.e., 
a~Pexp[—a/a^) with a cross-section weighted average ra- 
dius (summed over both components) of 0.092 /xm. The 
relative numbers of each grain type are such that the dust 
completely consumes slightly "super-solar" abundances in 
carbon and silicate: 400 ppm for C/H and 40 ppm for 
Si/H. We refer the reader to Whitney & Wolff (2000) for 
a more detailed discussion of the dust model. 

With the larger grains the amplitude of the variability is 
unaffecitcxl at low inc;linations because the direct starlight 
is much brighter than the scattered light. For i = 82°, we 
found larger variability (AV ~ Imag, A7 ~ 0.5mag) than 
for the simulation using KMH grains due to the increased 
beaming of the scattered light by the forward throwing 
grains. 

Increasing the phase function asymmetry through larger 
grains increases the photometric variability, but it is still 
less than observed. Therefore, we have investigated other 
models which have larger spot/star luminosity ratios and 
different spot parameters. Figure 3 shows a simulation in 
which we obtain AV ~ 1.45mags and A/ ~ l.OSmags, 
in agreement with our observations. This simulation has 
T, = lO^K, = 3000K, 6*, = 20°, and = 60°. We can 
obtain the same amplitude of variability with a model that 
has = 3800K and = 2 x 10*K. As the central star in 
the HH30 system is not observed directly, its effective tem- 
perature is somewhat uncertain and may be cooler than 
the T-t = 3800K, MO spectral typo determined by Kcnyon 
et al. (1998). The large spot/star luminosity, required to 
match our observations, could result in considerable veiling 
allowing for a later spectral type (see Kenyon et al. 1998). 
The spot size and latitude may be constrained through de- 
tailed modeling of time series images of the scattered light 
disk (Stapelfeldt, et al. 1999a). 

4. SUMMARY 

We have presented VRI photometric observations of 
HH30 that show it is highly variable over timescales of 
a few days (AF ^ 1.4 mags, AiJ ~ 1.0 mags, AJ ~ 
1.1 mags), with hints of periodicity at 11.6 or 19.8 days. 
This contrasts with the recent work of Stapelfeldt et al. 
(1999b) whose HST observations suggest a characteristic 
timescale on the order of 15 years, possibly related to in- 
homogeneities in the disk. The much shorter timescale 
variability in our data indicates a stellar origin for the pho- 
tometric variability we observe. High resolution imaging is 



required to determine whether the variability we observe 
is related to the morphological variations in the scattered 
light images. 

Comparing our findings to classical T Tauri stars, 
HH30's variability is among the largest yet reported. We 
have modeled the variability in the context of the magneto- 
spheric accretion model with a single hot spot on the stel- 
lar surface. The large variability requires a large spot/star 
luminosity ratio and that the circumstellar grains are 
larger and more forward throwing than interstellar grains. 
Variability due to accretion hot spots is often stochastic 
(e.g., Herbst et al. 1994) and while our modeling has 
adopted a single stable hotspot, further observations are 
required to better sample the lightcurve and determine 
changes in the spot sizes and temperatures. Our model- 
ing of the circumstellar dust is consistent with many other 
studies that are providing evidence for grain growth in 
dense environments of T Tauri stars. 

While our observations and modeling support the mag- 
netic accretion hypothesis for HH30; alternative sources 
of the variability may be orbiting dust clouds, inhomo- 
geneities in the disk structure close to the star, a binary 
star, and stellar flares. Further monitoring of HH30 with 
sampling on the order of days will yield further insight 
into this system. If the variability is due to magnetic ac- 
cretion there will be specific time sequence variations of 
the phometry, polarimetry, and morphology of the scat- 
tered light disk due to illumination by hot spots. Even 
if the accretion is stochastic there will be specific correla- 
tions between the photometry and polarimetry (Stassun & 
Wood 1999, Fig. 7). Rotational modulation will manifest 
itself in high resolution imaging studies through asymmet- 
ric brightening of scattered light images (Wood & Whitney 
1998; Stapelfeldt et al. 1999a). Lower resolution imaging 
should still be able to detect the asymmetric brightening 
via shifting of the photometric centroid of the system as 
pointed out by Stapelfeldt et al. (1999a). 
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Fig. 1. — HH30 VRI photometric variability. 



Table 1 
Parameters for Dust Grains 



Wave Band 


K (cm^/g) 


a 


9 


P 


y(KMH) . . 


220 


0.54 


0.44 


0.43 


/(KMH) . . . 


105 


0.49 


0.29 


0.70 


y (Large) . . 


174 


0.47 


0.61 


0.38 


/(Large) . . . 


134 


0.49 


0.58 


0.34 
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Fig. 2. — Single spot model with KMH grains, 'I'^ = lO^K, 'I\ — 3800K, 9s — 20°, (ps — 65°. Upper panel shows the variability for viewing angles of 
i — 30° and i — QO ^ lower panel shows the variability for the more edge-on viewing angle 2 — 82 appropriate to HH30. Solid line is the Y simulation 
and the dashed line is the / simulation. 
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